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Abstract 

Savannas are dynamical systems where grasses and trees can either dominate or coexist. Fires 
are known to be central in the functioning of the savanna biome though their characteristics are 
expected to vary along the rainfall gradients as observed in Sub-Saharan Africa. In this paper, 
we model the tree-grass dynamics using impulsive differential equations that consider fires as 
discrete events. This framework allows us to carry out a comprehensive qualitative mathematical 
analysis that revealed more diverse possible outcomes than the analogous continuous model. We 
investigated local and global properties of the equilibria and show that various states exist for the 
physiognomy of vegetation. Though several abrupt shifts between vegetation states appeared 
determined by fire periodicity, we showed that direct shading of grasses by trees is also an 
influential process embodied in the model by a competition parameter leading to bifurcations. 
Relying on a suitable nonstandard finite difference scheme, we carried out numerical simulations 
in reference to three main climatic zones as observable in Central Africa. 

Keywords: Savanna; tree/grass competition; ecological gradients ; fires; periodic solutions; stabil¬ 
ity; impulsive differential equations (IDE); bifurcation; nonstandard finite difference scheme. 

MSC Classification: Primary 30A37, 92D40. Secondary: 37M05 


1 Introduction 

In savannas, trees and grasses typically coexist jl]. Fire is recognized as playing a major part in the 
dynamics of this biome. The nature of grass-tree interactions and fire regimes strongly vary along 
environmental gradients in tropical savannas. Fire is more intense in wet than in arid savannas, 
where lower water availability leads to lower grass, i.e. fuel load, production. Thus fire is expected 
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to control Tree-Grass dynamics in wet savannas |2]. But two hypotheses for Tree-Grass coexistence 
have been introduced during these last decades. First, Walter (1971) j3] proposed the idea that 
trees and grass exploit two different rooting niches. Grasses are rooted in superficial soil layers and 
first use the incoming water, whereas tree roots are situated in subsoil, so that trees could grow only 
where enough water reached deeper soil horizons. This idea was developed analytically by Walker 
and Noy-Meir (1982) |3] using a Lotka-Volterra theory of co-existence between competitors. The 
second hypothesis says that grass-tree coexistence is driven by limited opportunities for seedling to 
escape both droughts and flame zone into the adult stage (Hochberg et al. 1994 |^; Higgins et al. 
(2000) [6]). In areas where tree seedlings succeed to establish in spite of competition with grasses, 
they are burnt by frequent grass fires (Higgins et al. 2000 |^). 

Savannas fires are frequent, up to occurring every 1-5 years in wet savannas (Frost & Robertson 
1985) though the fire return time is usually a decreasing function of mean annual precipitation. Fuel 
load made of dead aerial grass parts typically ranges between 2 and 10 t.ha“^ of dry matter (DM) 
(Lacey et al. 1982 [7]; Stronach & Mac-Naughton 1989 |8]; Menaut et al. 1991 |9]; Mordelet 1993 
|10| ) and flame height is usually 2-3 metres high (Frost & Robertson 1987 |11)L Although the fire 
burns most or all the aboveground grass biomass, the large underground root systems of perennial 
grass species enable most of the tufts to survive even the most intense fires and to rapidly establish 
new shoots before the onset of the rainy season. In contrast to grasses, trees which are less than 
2m height may either succumb to fire or have to resprout from roots and have their growth delayed 
(Bond and Midgley, 2001 [12] ). Mature trees (> 8m) and shrubs beyond 2m are more fire resistant 
and only experience partial die-back (Menaut &: Gesar 1979 |13j : Gillon 1983 |14|1. Early fires (in 
the beginning of the dry season) are less violent than late fires and have a lower impact on tree 
regeneration (Abbadie et al. 2006 |15)L 

Africa is a land of extreme contrasts in rainfall distribution and the time of year during which 
rainfall occurs (Janowiak 1987 |16|L When soil resource supply is temporally variable, trees and 
grasses will experience two distinct phases of resource availability: pulse periods when resources 
are high and most growth and biomass accumulation (fuel load) occurs, and inter-pulse periods 
when resources are too low for most tree and grass to take up and most mortality due to resource 
deficits takes place (Goldberg &: Novoplansky, 1997 |17| : Noy-Meir, 1973 |18]). Hence essential 
resource availability (e.g. water) is discontinuously available and the availability of these resources 
impact the ecosystem as discrete pulse events interspersed among long periods of limited resource 
availability (Schwinning et al. 2004 |19)1. 

Fires are sudden event that consume trees and grass biomass (Scheiter 2008 |20|). The broad 
objective of this study is to examine the influence of pulse events with regard to fires impact on the 
Tree-Grass dynamics along the rainfall gradient in Africa. Tree-grass savanna models can not be 
studied without the important role of fires (Tilman 1994 [21]; Higgins et al. 2000 [22]; Sankaran et 
al. 2004 [23], 2005 |24| : D’Odorico et al. 2006 [25]; Accatino et al. 2010 [26]; Beckage et al. 2011 
|27j : Staver et al. 2011 [28]; Yatat et al. 2014 |2^ and Tchuinte et al. 2014 |30j). This paper extends 
our earlier work (Tchuinte et al. 2014 |3^) where we consider a continuous tree-grass interaction 
model that featured a fairly generic family of non-linear functions of grass biomass to model fire 
intensity and its impact on tree. We have shown that the continuous model is able to predict a 
variety of dynamical outcomes. Notably, the number of equilibria featuring Tree-Grass coexistence 
depends on the characteristics of fairly generic Monod functions used to model the fire impact on 
tree dynamics. Moreover, we have shown that various bistability situations occur among forest, 
grassland and Tree-Grass (i.e. savanna) equilibria (for more detail see Tchuinte et al. 2014 |30)1. 
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Of course, in practice, fires are not continuous. Recent studies of the interactions between fire and 
vegetation are based on stochastic approaches because of the random and unpredictable nature of 
hre occurrences (D’Odorico et ah, 2006 [23; Beckage et al. 2011 |27]b 

In section 2 we will present the model with pulse hres. The theoretical analysis is developed 
in section 3. We show that the system admits four equilibria among which two trivial equilibria 
(the bare soil and the forest equilibria), and two periodic equilibria (the periodic grassland and 
the periodic savanna equilibria). We show that there are various bistabilities: between forest and 
grassland; between forest and savanna. Local and global stabilities are distinguished using classical 
tools such as Floquet multipliers and comparison theorem. We highlight thresholds that summarize 
the dynamics of the model and explain the theoretical meaning of these thresholds. Prior to illustrate 
our theoretical results numerically, in section 4, based on the scheme developed in |29) . we develop 
a reliable nonstandard hnite difference method (NSFD) that preserves the qualitative properties of 
the system (Anguelov et al. 2012 |31) . 2013 |32| . 2014 [33] )• Section 6 concludes the paper. Some 
mathematical details are included in appendices. 


2 The mathematical model 

In savanna environment, hre intensity is tightly linked to dried grass biomass that remains during 
the dry season (Higgins et al. 2008 )22)L During the last decades the effects of hre on vegetation 
dynamics have been studied (Scholes and Walker 1993 |34j : Higgins et al. 2000 |3)- Most of the 
models associated to or derived from these studies are ordinary differential equations (ODE) which 
assume that hres occur continuously with a hxed frequency. However hres are sudden event that 
consume grass biomasses and kill or harm tree seedlings (Scheiter 2008 |20)). The season of burning 
and the time between recurring hres determine trees and grass physioniomies in most ecosystems 
and especially in the savanna biome (Thonicke et ah, (2001) |35|L In this paper we present a 
new Tree-Grass model that aim to contribute to our understanding about how pulse hre shapes 
vegetation dynamics in hre-prone savanna-like ecosystems. We consider hre as discrete events and 
derive the following impulsive differential system 

’ f=-'co(l-^)-Sc»G-,raTG.' 

^ it ^ in 

< AG{tn) = G{t+) - Gitn) = -XfcGitn), 

AT{tn) = T{tt) - T{tn) = -XfMXfGG{tn))T{tn), 

G{4) = Go, 

. T{4) = To, 

where, 

• T and G are tree and grass biomasses respectively. 


,t — tn,Tl — 1,2,..., N-j- 


( 1 ) 
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• r = j is the period of time between two consecutive fires, and / is the frequency of fire, 

• Nr is a countable number of fire occurrence, 

• tn = nr, n = 1, 2,..., Nr are called moments of impulsive effects of fire, and satisfy 0 < < 

t2 < ... <tk< tNr ) 

• U}{G) is a generic non-linear functional which expresses fire intensity as an increasing function 

of grass biomass. Other than smoothness, it satisfies the following three conditions: (i) fires 

spread if and only if fuel is available (w(0) = 0), (n) fire-impact increases with fuel available 

(w(G') > 0, w (G) > 0), and {iii) there is boundary effects lim uj{G) < 1. 

G—>-00 

Other parameters used are listed in the following table. 


Table 1: Parameter symbols and names used to initialize the model 


Symbol 

Parameter name 

Units 

IG 

Grass biomass production per unit of grass biomass per year 

yr-i 


Grass biomass loss by herbivory (grazing) or human action 

yr-i 

Kg 

Garrying capacity of grass biomass 

t.ha“^ 

o 

II 

Additional death due to grass-grass competition 

ha.t“^.yr“^ 

^fG 

loss of grass biomass due to fire 

- 

It 

Tree biomass production per unit of tree biomass per year 

yr-i 

S'j' 

Tree biomass loss by herbivory (browsing) or human action 

yr-1 

Kx 

Garrying capacity of tree biomass 

t.ha“^ 

7t 

Additional death due to tree-tree competition 

ha.t“^.yr“^ 

XfT 

loss of tree biomass due to fire 

- 

ITG 

grass mortality due to tree/grass competition 

ha.t“^.yr“^ 


We suppose that solutions of Q is right continuous at tn, n = 1,2,Nr, that is G{t^) = 

lim G{tn + h) = Gitn) and Titn) = lim T(tn + h) = Titn), where G(t^) and Titn) are the 
h^0+ h^0+ 

biomass values for grasses and trees instantly after impulsive fire. Immediately following each fire 
pulse, system Q evolves from its new initial state without being further affected by the fire scheme 
until the next pulse is applied. In agreement with empirical experience (Abbadie et ah, 2006 jlSjh 
we assume that the level of destruction of the tree biomass depends on the available grass biomass 
through a;(G). 


3 Theoretical analysis 

System Q belongs to basic theory of IDE (Bainov 1993 |36 ) ) and their applications in Ecology. 
IDEs generally describe phenomena which are subject to abrupt or instantaneous changes. Model 


4 








Q derives from the family of impulsive Kolmogorov-type population dynamics in the theory of 
mathematical biology which the general form is given by 


dx 

dt 


= Xi{t)fi{t,x{t)), t / 4 , 


( 2 ) 


[ Xi{t'l) = Iik{tk,x{tk)), t = tk 

where Xi{t) represents the density or size of species Xi at the time t, x{t) = {xi{t), X2{t) , ...,Xn{t)), 
f{t,x) is an n-dimensional real functional defined by 

= {fi{t,x),f2{t,x), ...,fn{t,x)) 

and Ikitk, x) is also an n-dimensional real functional defined by 


Ik{tkjX') — (/ifc(4) ^))-^2fc(4) 3^)) ■ • (4) 3 ;)). 

This family of models has recently attracted the attention of several authors (Ezi, m, eh], m, 
|4l) . |42j . )43| ■ and the references cited therein). The main study subjects are the permanence, per¬ 
sistence and extinction of species, the local and global asymptotic stability of systems, the existence 
and uniqueness of positive periodic solution and almost periodic solution, and the bifurcation and 
dynamical complexity, etc. However, in all models investigated in the literature, authors generally 
do not consider the non-linear impulses i.e. non-linear form of the function Ik{tk,x). They mostly 
focused on the quasi-linear impulses. Since for the continuous model (see Tchuinte et al. 2014 |30|) 
the non-linear shape brings a wealth of possibilities, notably for the existence of various positive 
Tree-Grass equilibria. We retain this option although it could render the model difficult to study. In 
our model, we consider a generic non-linear functional response U}{G), which expresses the causality 
between grass biomass and fire intensity as to model the impact of fire on the woody biomass. 
Depending of the fuel accumulation (grass biomass i.e. G), U}{G) could take the general sigmoidal 
G" 

a non-linear increasing function. In this regard, since our Tree-Grass model does not contain more 
than two populations in competition, it may appear to be simple mathematically at first sight, but 
it is, in fact very challenging and complicated due to the nonlinear impulsive functions. 

Prior to analyzing the above model Q, it is important to show positivity and boundedness for 
solutions as they represent biomasses. Positivity implies that the populations survive and bound¬ 
edness may be interpreted as a natural restriction to growth as a consequence of limited resources. 
Then, model a requires that trajectories remain positive and that trajectories do not tend to 

infinity with increasing time. 

xq rj' 

Set Xq = — and Yt = —, where tq = 'yc ~ 4go and tt = 'yT — 4r are net primary production 
d'G d-T 

of grass and tree biomasses respectively. The following lemma holds. 

Lemma 3.1 When 'yc > 4go CLud 'yx > St, the eompaet 

B = S^{G,Tf G Kl/G < Xg = Kg (1 - ,T < Tt = (1 - I 

is positively invariant and attracting for the system 0 - Note that solutions of system 0 are 
bounded below by 0, and above by the carrying capacities of grass and tree biomasses. 


0 > 0, or some equivalent form. However, in our study uj{G) is principally treated as 
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Proof: See appendix A, page 

When t 7 ^ nr, the right-hand side of system Q is locally Lipschitz continuous on B. Thus, 
system Q has a unique solution. 


3.1 Equilibria 

System Q has constant and periodic equilibria. 

First of all, it is obvious that Eqq = (0;0) and Eqi = (0, It) are "trivial" equilibria of system Q. 
While Floo represents the bare soil, Eqi is the constant forest equilibrium, like in Tchuinte et al. 
2014 [30]. 

Now, let us show the existence of periodic solutions of the impulsive system 0. The existence 
of the periodic grassland equilibrium depends on the following threshold 




Ge 

0,pulse 



T 


rc 

(1 - -^/G 


The following theorem holds. 


Theorem 3.1 

jG. 


(Semi-trivial periodic equilibrium) 
when > 1, System (1) has a periodic grassland equilibrium Eiq = (Ge(t);0), where 


Geit) 


Xg 


{(1 - \fG)E'G'^ - l}e^G{t-nT) 


t G [nr, (n -|- l)r[, n = 0,1, 2,... 


Proof: See appendix B, page[^ 

Let us now show that there exists an unique positive periodic Tree-Grass equilibrium {G* {t)]T*[t)). 
We set 


m{t, nr, T^nr)) 


rG{t — nr) lu 

hT 


1 + 


r(nr) 

Yt 


1 

gr-rO-nr) _ ^ 


X(t,nr,r(nr)) 


_ mit^nr^Tinr)) 

^ 1 


(3) 

(4) 


and. 


where 


G* = 


[(1 - \fG)x{{n + 1)a nr, T*) - 1] 


tJ^G Xiu, nr, T*)du 

yT{(l-A/rw(A/GG*))e^^--l} 


, n = 0,l,2,.... 


_ 

7?* 


- 1 ) 

rr 


-In 


1 


Ge(r) = Ag 


^ V ^ “ A/Ta;(A/GG'e(r)) ^ 

(1-A/G)e’'«"-1 


e'-GT _ 1 


(5) 

( 6 ) 
(7) 
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Theorem 3.2 (Uniqueness of the non-trivial periodic equilibrium) 

When > 1 md "T^opuZse ^ then system (1) has a unique positive Tree-Grass periodic 

equilibrium = {G*{t)]T*{t)), where 


G*{t) = 


x{t,nT, T*)G* 


and T*{t) = 


^rT{t-nT)rp* 


1 + x{u, nr, T*)du 

nr < t < (n + 1 )t, n = 0,1, 2 ,such that G*{t) = G* and T*(r) = T*. 


1 + - 1 ) 

Yt 


Proof: See appendix C, page 31 


3.2 Local stability of the "trivial" equilibria 

Set 

noi = ^ X and iZQ,noi = pulse ■ 

rr JTG ^ V '^'01 / 

Theorem 3.3 (Local stability of constant equilibria) 

1. Eqq = (0; 0) is always unstable. 

2. IfTZoi < 1, then the forest equilibrium Eqi is locally asymptotically stable (LAS) (similarly as 
in the continuous model (see Tchuinte et al. (2014) }S0^ f 

3. If TZqi > 1 and IZq^tioi < 1? then Eqi is LAS. This situation is specific for the impulse model. 
The continuous model does not imply the stability of the forest when IZqi > 1. 

4 . IfEoi > 1 and 'T^o.'R.qi > then Eqi is unstable. 


Proof: See appendix D, page 37 


3.3 Local stability of periodic equilibria 

We begin to investigate the local asymptotic stability of the periodic grassland equilibrium of system 
Q. To complete this subsection, we show the local stability of the periodic savanna equilibrium. 
Recall that 


X 


Ge 

0,pulse 


Tg 


-In 

T 


1 - A 


fG 


and IT, 


rT 


O^pulse 


-In 

r 


1 — XfTUj{XfGGe{T)) 


where 


Ge{t) = Xg 


{(l-A/G)e^GT_l}erG(t-nr) 

{(1 - A/G)e'’G^ - l}e^G(i-nr) Xj^e'^G-r 


t £ [nr, (n + l)r[, n = 0,1, 2,.... 


~ 

Theorem 3.4 If > 1 and ^ periodic grassland equilibrium Eiq = 

(Ge(t),0) is locally asymptotically stable. 


7 
















Proof: See appendix E, page 39 

Now, we investigate local properties of the periodic savanna equilibrium. Set 




TT 


0,stable 


and 


1^^/ 1 

Vl “ ■^/tw(A/gG*) 

1 


- — 
'^0,7?.oi -Tp 
/■Cl 


1 


1 - 


1 


01 


n, 


{),stable 


stable 


- ^^7^01 + ^Ge 


K 


O^pulse 




where G* and T* are defined in (j5| and (j6j) respectively. The threshold stable represents the 
net production of tree biomass relative to nre-induced tree biomass loss at the mixed Tree-Grass 
equilibrium. 

The following theorem holds. 

Theorem 3.5 When TZ*^ > 1 and TZ**^ > 1, the savanna periodie equilibrium = 
{G*{t)]T*{t)) is LAS. with 


G*{t) = 


x{t,nT, T*)G* 


and T*{t) = 


^rT{t-nT)rp* 


1 p,gG* x(u, nr, T*)du 
nr <t < {n + 1)t, n = 0,1, 2,..., where the expression of x is given in 0- 


1 + —(e’-Vt-nr) _ j) 

Yt 


Proof: See appendix F, page 41 


The local stability is sufficient when there are multiple stable states. However, for an unique 
equilibrium, the global stability is necessary to ensure that all trajectories converge to the equilib¬ 


rium. 


Remark 3.1 We compute G* and T* using ^ and ^ and a specific command ("fzero") in matlab 
which determines the fixed point. This allows us to obtain ’^os/afe/e- numerical approxima¬ 

tions of G*{t) to have 1Zl*stahie- 


3.4 Global stability of equilibria 

In this section, we investigate the global stability of the forest equilibrium and the periodic grassland 
equilibrium. The following theorem holds. 

Theorem 3.6 (Forest equilibrium GAS) 

The Forest equilibrium Eqi = (0; Yr) is globally asymptotically stable when < 1. 

Proof: See appendix G, page|44| 

The global stability of the periodic grassland equilibrium is given in the following theorem. 












Theorem 3.7 (Grassland periodic equilibrium GAS) 

If^o^yaise ^ ^ '^ 0 pulse ^ then the grassland periodic equilibrium ^lo = {Geit)]0) is globally 

asymptotically stable, where 


Ge{t) = Xg 

n] 


{{1 - XfG)e^^^ - l}e^G(t-nT) 

{(1 - Xf^e^G-r 

Ge _ re 


0,pulse 

— In 


1 


T \l-XfG 
Proof: See appendix H. 


and TZ, 


t G [nr, (n + l)r[, n = 0,1, 2, 
rr 


0,pulse 


-In 


1 


^ yl “ "^/Tw(A/GGe(r))y 
Remark 3.2 Gonsider G*{t) < Ge(r). We show that stable- Therefore 


• '^0,stable < 1 '^0,pulse < 

• '^0,pulse > 1 TIq, stable > 


Thresholds and their ecological meaning are recalled in the following table 


Table 2: The thresholds and their ecological meaning 


Thresholds 

Ecological meaning 

K„, = ^2 X ir 

rr ITG 

the net primary production of grasses relative to the grass 
production loss due to tree/grass competition 

throughout their life at the close forest equilibrium 

^Ge _ ec 

'^{),pulse 1 / 1 \ 

Vl - A/g j 

It is the net primary production of grasses after fire 

T^O.noi = X^fpulse (^1 “ 

It is the mixed threshold of the two previous thresholds 

TZ* — 

‘^f),pulse ^ ^ \ 

r y 1 - XfTUj{XfGGe{T )) j 

represents the net production of tree biomass relative 
to fire-induced biomass loss at the period of fire 

at the grassland equilibrium 

TZ* — 

'^0,stable ^ ^ \ 

T y 1 - A/tw(A/gG*)J 

It is the net production of tree biomass relative 
to fire-induced biomass loss at the period of fire 

at the mixed Tree-Grass equilibrium 


We summarize all local and global properties of the impulsive model Q in table 
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Table 3: Long term behaviour of model Q 


Thresholds 

Equilibria 

Stable 

Unstable 

Case 

Ti-oi 

-rjGc 
^0,-pulse 

72.0,TZoi 

7? * 

'^0,pulse 

^0,sta6Ze 

stable 


<1 


- 

- 

- 

Eoo, Eoi 

Eoi (GAS) 

Eoo 

I 




<1 

- 

- 

Eoo, Eoi 

Eoi (LAS) 

Eoo 

II 



- 




Eio 

Eio (LAS) 




>1 




>1 

Eoo, Eoi 

Eoi (LAS) 

Eoo 

III 

< 1 


>1 

>1 


Eio, Eh 

Ell (LAS) 

Eio 






<1 

Eoo, Eoi 


Eoo, Eio 

IV 





<1 

- 

.Bio, Eh 

Eoi (GAS) 

Ell 








Eoo, Eoi 


Eoo 






- 

- 

Eio 

Eoi (GAS) 

Eoi 

V 


<1 


- 

- 

- 

Eoo, Eoi 

Eoi (GAS) 

Eoo 

VI 




<1 

- 

- 

Eoo, Eoi 

Eoi (LAS) 

O 

o 

VII 



< 1 




Eio 

Eio (LAS) 




>1 




>1 

Eoo, Eoi 

Eoi (LAS) 

o 

o 

VIII 



>1 

>1 


Eio, Eh 

Ell (LAS) 

Eio 






<1 

Eoo, Eoi 


o 

o 

o 

IX 





<1 

- 

Eio, Eh 

Eoi (GAS) 

Ell 






- 

- 

Eoo, Eio 


Eoo 

X 

> 1 






Eoi 

Eoi (GAS) 

Eio 


<1 


- 

- 

- 

Eoo, Eoi 

Eoi (GAS) 

Eoo 

XI 




<1 

- 

- 

Eoo, Eoi 


Eoo, Eoi 

XII 



>1 




Eio 

Eio (GAS) 






>1 


>1 

Eoo, Eoi 


Eoo, Eoi 

XIII 


>1 



>1 


Eio, Eh 

Ell (GAS) 

Eio 







<1 

Eoo, Eoi 


Eoo, Eio 

XIV 





<1 

- 

Eio, Eh 

Eoi (GAS) 

Eh 






- 

- 

Eoo, Eio 


Eoo 








Eoi 

Eio (GAS) 

Eoi 

XV 


Table 1^ gives all possible configurations for the impulsive system. According to the values taken 
by the thresholds, it may be possible to anticipate the long term behaviour of the system. Since many 
configurations are possible, it is essential to highlights the parameters that may have an important 
impact on the thresholds. In the next sections, we briefly present the numerical algorithm we have 
chosen to perform numerical simulations and present some results emphasizing the importance of 
the competition parameter 7 tg- 
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4 The numerical algorithm 


In the previous section, solutions were searched in form of analytical expression. However, many 
impulsive differential equations can not be solved in this way or their solving is more complicated 
in the mathematical point of view. 

A nonstandard numerical scheme for solving the impulsive differential equation is built. The 
nonstandard approach relies on the following important rules: the standard denominator At in 
each discrete derivative is replaced by a time-step function 0 < ip{At) < 1 ; such that (p{At) = 
At + 0{At)] the nonlinear terms are approximated in a non local way; for instance the nonlinear 
term T{tn)G{tn) in the problem can be approximated by For an overview and some 

applications in Biology of the nonstandard finite difference method see for instance (Anguelov et 
al., 2012 |31| : Anguelov et ah, 2014 |33ji. 

The nonstandard approximations for system Q are given by 


Qn+l _ QT 

(pi{At) 

2^72+1 _ 

ip2{At) 


(7G - Sgo)G^ - 


(7t - St)T^ - 


( 8 ) 


and 


where 


and 



= (1-A/g)G-+i, 

= {l-XfTCo{XfGGn)T^+\ 

g(7G-<5Go)Ai _ 


q(At) = 

72{At) = 


IG - <^G0 

g(7T—(5T)Ai _ ^ 
'yn^ — HT’ 


(9) 


( 10 ) 

( 11 ) 


Scheme ([^b) with the time-step function (11) is an exact scheme, between each fire event. Similarly, 
when ’’Ytg = 0) the scheme (|^a) with the time-step function (10) is also an exact scheme, between 
each fire event. Altogether the numerical algorithm (§-(§ is positively stable and elementary 
stable i.e. it preserves equilibria and local properties of each equilibrium of system (Q. Thus at 
least locally, we are sure that schemes 0-([^ replicate the dynamics of system (Q. 


5 Numerical simulatious aud discussiou 

5.1 Environmental setting 

To illustrate our analytical results and highlight important ecological parameters, we will perform 
some numerical simulations. In fact, from the ecological point of view, parameters can change 
drastically according to the environmental features. For instance, in Cameroon, three different zones 
can be particularly highlighted. The first zone is Region 1 (Rl) where the biomass production is low 
(the Mean Annual Precipitations (MAP) is less than 650 mm by year). The second zone is Region 2 
(R2) (650-1100 mm/yr). It has more biomass production by year. The last zone (R3) (1100 — 1800 
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mm by year) where savannas are observed in some cases in the immediate vicinity of forests (Favier 
et al. 2012 )44j). Model 0 has 11 parameters. Eight of them (yc, iFc, <5 go, a, Tr,-fC-r, J-r, 7 ^) 
specify vegetation growth while the others are related to the fire characteristics. Parameter values 
used are based on literature sources. 



Figure 1: Overall vegetation map of Cameroon from Letouzey (1985). The three regions Rl, R2 and R3 are 
indicated by black rectangles. 
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Table 4: Vegetational Data 


Parameters 

Semi — arid 

Rl 

Site name 

Mesic 

R2 

Humid 

R3 

Units 

Sources 

Kg 

2-5 

8-10 

10-20 

t.ha“^ 

IISIITS] 

7g 

0.4i*i - 1.5i**^ 

1.5-3 

3 - 4.6^***) 

yr-i 

i*/ gni, [HI 

5go 

0 

0 

0 - 0.9^^*> 


(4*) [47| 

^fG 

0.1-0.9 

0.1-0.9 

0.1-0.9 

- 

assumed, see also |15| 

a 

2 

2 

2 

t.ha“^ 

Assumed 

e 

2 

2 

2 

t.ha“^ 

Assumed 

Kt 

10-25 

25-60 

60-115 

t.ha~^ 

m 

7r 

O 

1 

o 

0.9- 1.2 

1.2 - 7.2^*’) 

yr“^ 

W jig, [49] 

St 

0 

0 

0-0.015^®*) 

yr"^ 

(5*) [5] 

XfT 

0.1-0.5 

0.1-0.5 

0.1-0.4 

- 

Assumed 

T 

> 10 

2-8 

0.5-2 

yr 

Overall expert-based knowledge, 

in addition to |13[ I44| 

'yTG 

(-0.01) -0.03 

0.01-0.08 

0.03 - 0.09 

ha.f^.yr”’^ 

ESI, ESI 

(after reinterpretation) 


Note 5.1 Range values for parameters (Kq, kr, ^fG, ^fT, r, 'yra) used for the simulation runs. Early fire destroy 
only 25% (Abbadie et al, 2006 \15^ ) of grass biomass, while late fires can destroy up to 90% of biomass. 


Our numerical analysis focuses on Cameroon as part of Central Africa where we find a summary 
of African natural conditions, from humid equatorial climate near the Atlantic Ocean, up to the 
arid Sahelian tropical climate in the region of Lake Chad. 

The first site is Rl. It corresponds to semi-arid zones. Grasses may be dominant while trees are 
generally of low stature. Trees are resource-limited, and the resource competition with grasses and 
between trees are the key factor determining savanna existence (Baudena et ah, 2014 [50]; Tchuinte 
et al., 2014 |30|). In Cameroon, Rl (small black rectangle in Figure]^ corresponds to the dry 
tropical climate of the extreme North Country, from Kaele in Maroua and Mora, and from Yagoua 
to Kousseri, Makary and Lake Chad. At the edges of Lake Chad, there is only 3 months of rain 
with 500 mm/yr. 

The second region, R2, represents a mesic zone. Here the MAP varies from 650 mm/yr to 
1100 mm/yr. In Cameroon it is located at the Northern part, from Adamawa to the Mandara 
Mountains (see the middle black rectangle in fig [^. There are two seasons covering the whole 
of Adamawa plateau, from Banyo to Ngaroundere and Meiganga. It has been argued, that mesic 
savannas are unstable compared to forest and disturbance-dependent with respect to fires (Sankaran 
et al., 2005 jl^), which prevent tree invasion, because they occur regularly during the dry season. 
Grasses benefit from fire because they recover faster than trees after fires, and profit of open spaces 
to growth. Thus grass-fire feedback is a characteristic feature that leads to savanna or grassland 
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persistence. 

High rainfall occurring in wet African savannas directly reduce the role of water as limiting 
factor. The last region (R3) that we are interested in corresponds to humid areas where MAP 
varies between 1100 mm/yr and 1800 mm/yr. In this region, the high water resource available 
enables high fuel (grass biomass) production and therefore fires are more frequent and of greater 
impact on seedlings. The grass-fire feedback in R3 leads to a bistability of savanna and forest, as 
shown using a simple continuous models (e.g. Tchuinte et ah, 2014 m; St aver and Levin, 2012 
|51| 1 and evidenced from remote sensing data by Favier et al. (2012) |44) . In Cameroon, R3 (see 
the big black rectangle in Figure Q encompasses two sub-zones: a sub-zone of transition between 
equatorial and tropical climates, and a sub-zone which corresponds to the equatorial climate itself. 
Concerning the first one, the MAP is between 1100 mm/yr and 1500 mm/yr. It is observed from 
Bafia to Bertoua, Batouri and from Yoko to Betare Oya, Garoua Boulai'. The second sub-zone 
is a site covering the entire South of the Country from Yaounde (1564 mm/yr) to Yokadouma, 
from Ebolowa to Ambam, Mouloundou and Ouesso (Congo). It extented near the Gabonese borber 
(1700mm). In both the two sub-zones there are four distinct seasons (two dry seasons alternating 
with two wet seasons with unequal intensity). At the South Cameroon near the Gabonese border 
(11 months of rainy season), there is a close canopy forest (see hgQ. Above R3 (at the highest 
end of the rainfall range which is 2000 mm/yr), fires are totally suppressed and only forests are 
observed, since grass growth is inhibited by tree shade. 

In the next section we present some numerical simulations. The fundamental tasks in studying 
disturbance are to discriminate between fluctuations that are extraordinary and those that are usual 
(McNaughton 1992 |52)1. In our simulations, we assume that the ecological system is not impacted 
by Human and Animals (grazing and browsing), i.e. 5go = 5t = 0- 

5.2 Simulations to illustrate bifurcations due to jtg in regions Rl, R2, and R3 
5.2.1 Simulations in region Rl 

According to table we choose the following values of parameters for region Rl: 

Table 5: Parameters values related to figure 2 


Kg 

7G 

<5go 

Kx 

It 

5x ot 

T A fT 

A/g 

4 

0.7 

0 

14 

0.75 

0 2 

12 0.9 

0.5 


Taking 'jtg = [—0.01,0.01,0.03,0.051], and, using Table we obtain Table Figure illus¬ 
trates also the expected behaviours. 
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Table 6: Thresholds Table related to Table and Figure 


Panel 

’^01 

■jyGs 
^0,pulse 

T^o,noi 

V* 

^0,pulse 

7?* 

^O^stable 

^O^stahle 

Case 

a 

- 

> 1 

> 1 

> 1 

> 1 

> 1 

XIII 

b,c 

> 1 

> 1 

> 1 

> 1 

> 1 

> 1 

XIII 

d 

< 1 

> 1 

- 

> 1 

- 

- 

V 


Rib: yT£;=-«-«l 



Rib: 



Rib: yT£.=0.01 




Figure 2: Phase diagrams in Rl. Panels a and b (c and d) show the impact of the tree/grass competition parameter 
(see panels panel c and d). When 7 tg increases and exceeds a critical value, the attracting state of the system shifts 
from periodic savanna to a forest equilibrium (in fact, dense thickets). 


In semi-arid areas, Tree-Grass interactions are predominantly influenced by competition for soil 
water (e.g. Walker et al. 1981 [53]). However, shading by tree foliage under arid climate can also 
increase grass production under the tree crown (Abbadie et al. 2006 |15) i and more generally can 
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increase the water budget below the canopy (Barbier et al. 2008 |54|1. Hence the influence of trees 
on grasses can range from facilitation to competition. Figure shows the influence of tree/grass 
interactions in Rl. For higher values of 7 tG) trajectories converge to the forest equilibrium (see 
panel d) corresponding, in fact to dense thickets. When yrc is small, and even negative (positive 
effect on the grass biomass), the system converges to a periodic Tree-Grass coexistence equilibrium 
with fairly large amplitudes in the tree biomass (compare panels a and b), or with small grass 
biomass (see panel c). Here 'yxG is an influential parameter since it permits a transition from 
savanna to forest. This result joins those of Sankaran et al. (2005) |23] which argued that in Rl, 
savannas are stable in the sense that tree biomass and cover are primarily limited by resources (see 
panels a and c in Figure]^. Therefore, the competition parameter ■'yj'G is an important driver of 
Tree-Grass dynamics in Rl where fires return time are typically higher than 10 years and therefore 
are not necessary for grass-tree coexistence. Our Tree-Grass impulsional model shows that in Rl, 
there is only a stable periodic grass-tree equilibrium or a stable forest equilibrium (see panels a and 
d in Figure 1^. 


5.2.2 Simulations in region R2 

Let us consider the following values of parameters according to table 


Table 7: 


Parameters values related to figure 


Kg 

7G 

<5go 

Kj' 

7t 

5t 

Oi T \ 

A/g 

8 

1.9 

0 

30 

0.9 

0 

2 5 0.5 

0.6 


According to Table and Table with ''^tg = [0.01,0.02,0.03,0.055], we derive, in Table the 
behaviours of the Tree-Grass system. See also Figure page |17| 

Table 8: Thresholds Table related to Table and Figure 


Panel 

77-01 

-^Ge 
'^O,pulse 

77-0,7eoi 

'7?* 'ip** 

^0,pulse ^0,stable ^0,stable 

Case 

a,b,c 

> 1 

> 1 

> 1 

>1 >1 >1 

XIII 

d 

> 1 

> 1 

< 1 

>1 

X 


Figure l^below illustrates the bifurcation due to the competition parameter in Region (R2). The 
forest equilibrium is stable for ^tg = 0.055 a value in the upper range of plausible values (see panel 
d). When 'yxG decreases, the system converges to a savanna periodic equilibrium (see panel a, b, c). 
We note also that has an impact on the amplitude of the periodic savanna equilibrium and the 
maximal amount of grass biomass. 
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Grass biomass (t/ha) 


Figure 3: Phase diagrams in R2 with a period of fire of r = 5 years. 


Increasing the impact of fire on trees via an increase of the "A/t" coefficient, a bistability 
between forest and savanna occurs in R2 (see figure below). 

Table 9: Parameters values related to figure 


Kg 

7G 

<5go 

Kj' 

7t 

6 t 

ex T X j-j' 

A/g 

8 

1.5 

0 

30 

0.9 

0 

2 2.2 0.8 

0.5 


Using the same values for 7 tG) we derive Table 10 and Figure]^ 
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Table 10: Thresholds Table related to Table and Figure 


Panel 


■jyGs 
^0,pulse 

T^o,noi 

V* 

^0,pulse 

7?* 

^O^stable 

^O^stahle 

Case 

a, b 

> 1 

> 1 

> 1 

> 1 

> 1 

> 1 

XIII 

c >1 

> 1 

< 1 

> 1 

> 1 

> 1 

VIII 

d 

< 1 

> 1 

- 

> 1 

- 

- 

V 



Figure 4: Phase diagrams in R2 with a period of fire of r = 5 years (Table and a stronger impact 
of fire than in Figure 3 via XfT = 0.8. Panels a, b, c and d show the impact of the tree/grass 
competition parameter. 


Figure]^ shows interesting behaviours. In particular, in panel c, we obtain a bistable situation, 
where the system can either converge to the forest equilibrium or the periodic Tree-Grass equilib¬ 
rium, depending on the initial values. In that case, we don’t have analytic results that allow us to 
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know what are the basins of attraction of each equilibrium. That is why the use of a well fitted 
numerical scheme is of utmost importance, in order to capture this essential information. We will 
show other examples of bistability in the next section. 

5.2.3 Simulations in region R3 

According to table we first consider the following values for simulations in region R3: 

Table 11: Parameters values related to figure 5 


Kg 

7G 

5go 

Kx 

It 

5t 

ex T X j'j' 

A/g 

17 

4.5 

0 

45 

6 

0 

2 0.6 0.4 

0.4 


Using Table 15 with jtg = [0.03, 0.05, 0.07, 0.09] leads to Table 12 and figure )[^ 

and Figure 


Table 12: Thresholds Table related to Table 
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Panel 

77-01 

-^Ge 
''^0,pulse 

77-0,Teoi 

'ip* 'ip* '1^** 

^0,pulse ^0,stable ^0,stable 

Case 

a,b,c 

> 1 

> 1 

> 1 

>1 >1 >1 

XIII 

d 

> 1 

> 1 

< 1 

>1 

X 
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Figure 5: Phase diagrams in R3. This figure shows two equilibria: forest and periodic savanna. 


In the humid zone, the vegetation is intrinsically dominated by trees which exert competitive 
pressure on grasses, such that grasses are suppressed or even out-competed (Scholes and Walker 
1993 |34) . see panel d). Figure [^illustrates the effect of the competition parameter in R3. 

Note 5.2 Using realistic ranges for parameters, we show that in R3, a stable periodic savanna equi¬ 
librium may appear but also a stable forest equilibrium, for sufficiently high values of ^tg- However 
we cannot have a periodic grassland equilibrium. In R2, it is possible to have forest equilibrium and 
periodic savanna. 

Like in Region R2, bistable situations can occur in region R3. Let us first consider the following 
parameter values 
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Table 13: Parameters values related to figure 


Kg 

7G 

<^G0 

A/g 

Kx 

7t 

5t 

A/t 

r a jtg 

19 

3.1 

0.1 

0.5 

50 

1.5 

0.015 

0.6 

0.5 2 0.09 


Table 14: Threshold values related to figure 


’^01 '^o'pulse '^0,'R-oi T^*0,pulse '^0,stable '^0*stable ^ase 

< 1 > 1 - <1 - - II 


The mathematical analysis shows that there is a bistability between the forest equilibrium and 
a periodic grassland equilibrium (see line II in table [^. Panel a (pulse model) in figure [^illustrates 
two basins of attraction: one in favour of forest equilibrium; another in favour of the periodic 
grassland equilibrium (bistability as in line II in table [^. For the same values of parameters, the 
continuous model does not yield bistability (see panel b inj^: the only equilibrium is the forest 
equilibrium |30j . 


Bistability (pulse model) 



Grass biomass (t/ha) 


Monostability (continuous model) 



Figure 6: Comparison of the pulse vs. continuous models in reference to R3. The continuous model presents a 
forest equilibrium which is GAS (see panel b). In contrast, the discrete model shows two equilibria: the forest (as 
for the continuous model) and the periodic grassland (see panel a). Depending of the initial conditions the system 
converges to a stable periodic grassland or to a stable forest equilibrium. 


We have shown that the interspecific parameter, 'yj'G plays a great role in the dynamics. How¬ 
ever, fire period can also have an impact on the dynamics of the Tree-Grass system in R3 as in the 
two previous zones. Let us consider the following parameter values: 
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Table 15: Parameters values related to figure 


Kg 

7G 

<^G0 

A/g 

Kx 

7t 

5t 

A/t 

a jtg 

19 

3.1 

0.1 

0.5 

65 

1.5 

0.015 

0.6 

2 0.04 


According to Tables and 15 we derive, in Table 16, the different possible dynamics ot the 
Tree-Grass system: 


Table 16: Threshold values related to figure 


Panel 


■jyGs 
^0,pulse 


^{}^pulse ^O^stable ^O^stahle 

Case 

a 

> 1 

> 1 

< 1 

<1 

VII 

b 

> 1 

> 1 

< 1 

>1 >1 >1 

VIII 


It has been evidenced that in humid savannas, fire is necessary to establish the Tree-Grass 
coexistence equilibrium (Sankaran et al. 2005 |24)i. However, various bistabilities occur: between 
forest and grassland (see panel a in figure ; between forest and savanna (see panel b in figure . 
The system can shift from a to b (bifurcation) when the fire period increases. When r = 0.5 (two 
fires a year in a sub-equatorial context), we have a bistability case: with a forest equilibrium and a 
periodic grassland equilibrium (see panel a in figure which illustrates case VII in table [^. When 
the fire period increases by 10% i.e. r = 0.6, there is still bistability case, but the periodic grassland 
equilibrium is replaced by a periodic savanna equilibrium. This corresponds to line VIII in table 
and to panel b in figure 

In R3, both forest-grassland, and forest-savanna equilibria are predicted (see panels a and b in 
Figure]^. The fire period is a bifurcation parameter that shapes the Tree-Grass dynamics in R3. 


6 Conclusion 

Savannas are complex systems due to the interaction of trees and grasses which are frequently 
mediated by disturbances and notably by fires. The broad objective of this work was to develop a 
predictive understanding of Tree-Grass dynamics across rainfall gradients in Africa on the basis of 
a minimalistic model. This is done using specific features of three ecological contexts: semi-arid, 
mesic and humid. They represent different ecological conditions in terms of rainfall amount and 
deriving variations of most of the parameters used in the model. We formalize a new model of 
Tree-Grass interactions. The novelty of the model, with respect to other models (Staver et ah, 2011 
[28], Accatino et ah, 2010 |26] . |30| 1 is that fire is considered as discrete events with high or low 
return times. In addition, fire frequencies and fire intensities are decoupled. Discrete events are 
typically modelled by impulsive differential equations. We show here that this framework yields 
richer qualitative behaviours than continuous modelling (see figure]^. Several authors, in order to 
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BistabUity (t=0.5) 



Bifurcation (t=0.6) 



Grass biomass (t/ha) 


Figure 7: Fire mediated Tree-Grass mixtures and changes in their physiognomies according to the 
fire period. Panel a shows two bistable equilibria (a forest equilibrium and a periodic grassland 
equilibrium (see VII in table |^). From r = 0.5 to 0.6, the ecosystem changes. Panel b illustrates 
two bistable equilibria (the periodic savanna and the forest equilibrium (see VIII in table®). 


deal with the stochastic occurrence of fire have modelled fire through purely stochastic differential 
equations (D’Odoricco et al. 2006 |25) . Beckage et al. 2011 |27) . De Michele and Accatino 2014 
|55|). However results of those stochastic models are often obtained numerically by iterating the 
equations given parameter values and initial conditions and the authors have difficulties to verify 
mathematically that the qualitative properties (e.g existence of equilibria and their stability) of the 
model are preserved in their simulations. By contrast, at least qualitatively, our impulsive model is 
lend itself to a comprehensive qualitative analysis of the possible dynamic outcomes and properties 
of the system. For this raison, it is perhaps interesting to used impulsive framework to look at 
least qualitatively at the problem of predictability of discrete fire impacts in tropical savannas. The 
impulsive modelling of fire suggests ways for deriving from a minimal continuous fire model (e.g. 
Tchuinte et ah, 2014 |30) 1 more realistic discrete fire model. The theory of IDE allows for more 
detailed analyses of the system than stochastic differential equations. Further, using IDE technique, 
we can highlight mathematically thresholds that summarize the Tree-Grass interaction (see table 
and point out bifurcation parameters. 
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The impulsive model illustrates different kinds of dynamics which can be observed across of 
rainfall gradients in Africa. The model generates savanna equilibria in all the three regions. It 
means that trees are able to persist while not reaching 80-100% cover all over the wide range of 
rainfall considered. However, at high rainfall sites such like the boundaries of the tropical rainforests 
of central Africa, the vegetation is due to be dominated by trees and may even reach a closed tree 
cover in the absence of recurrent hres of low return times which are fostered by high grass production. 
On the other hand, trees can reach high individual biomass and exert intense light competition on 
grasses, such that grasses are suppressed or even out-competed (Scholes and Walker 1993 |34|1. This 
explains why forest may be stable in R3 and even in R2. However, tree cover can also facilitate 
the growth of grasses, i.e. 'Jtg < 0 in arid conditions (i.e. Rl) by limiting soil water losses from 
transpiration in the topsoil. In Rl, we show that the existence of savanna is mainly due to the 
competition parameter. Thus, only two equilibria are possible: a forest equilibrium and a periodic 
savanna equilibrium. The bifurcation from the forest to the periodic savanna is related to the 
resource competition parameter, 'yxG- This result joins those of Sankaran et ah, (2005) j24) that 
argued that in drier sites, savannas are stable in the sense that tree cover is intrinsically limited 
by resource (and this contributes to the low value of KT) and hre is not necessary for Tree-Grass 
coexistence. Recent modelling studies by Baudena et al. (2014) |50) conhrm that in semi-arid 
savannas, while trees are water-limited, the water competition with grasses is also a key factor 
determining savanna existence. 

We show that in all regions, the competition parameter plays a great role: it is a bifurcation 
parameter whatever the context. 

Grass-hre feedback principally occurs in mesic and humid areas. No savanna or grassland would 
emerge without this positive feedback in Regions R2 and R3. Using realistic parameter values, 
we show that in R2 (mesic area) forest and savanna are more present than grassland. In R3, in 
humid area, two bistability situations may occur: bistability between forest and grassland and 
bistability between forest and (periodic) savanna. In the equatorial climate in Southern Gameroon 
(Gentral Africa), there is forest-savanna contact in three sites inside R3: Akonolinga, Bertoua and 
Mbam-Kim (Youta 1998 [56]). Bistability of forest and grassland can be found in the equatorial 
and tropical climate of transitions (1100-1500 mm/yr) where hres are usually occurring every 1-5 
years (Frost and Robertson, 1987 |llj . Favier et al. 2012 jU). In this area and due to the large 
grass biomass, hame height is usually 2-3 m high (Frost and Robertson, 1987 |11]). Therefore, a 
severe hre could have a great impact on young trees/shrubs. In R3, water availability enables high 
fuel production. As a result, hre is severe and may occur as frequently as 0.5-1 yr. Thus, hre 
has a stabilizing role of grassland and savannas by preventing tree invasion on long time scales, 
freezing the forest-savanna boundary in a historical position (Gillon, 1983 |14jl. We show from our 
model that in R3 the Tree-Grass system can shift from bistability between forest and grassland to 
bistability between forest and savanna due to hre period. 

To conclude: (i) in all regions Tree-Grass competition is the most important parameter for 
Tree-Grass co-existence; (ii) in R2 and R3, hre period can also be a bifurcation parameter. Thus, 
the hre period (and hre intensities) and the competition parameters are the main determinant in 
the Tree-Grass dynamics; (iii) Modelling hre as pulse events provides more realistic situations than 
modelling hre as continuous events. 
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APPENDICES 

Some detail of the proofs of results associated with system 0 are provided. 


Appendix A: Proof of Lemma 3.1 


It is obvious that G = 0, and T = 0 are vertical and horizontal null-clines respectively. Then, no 
trajectory can cut these axes. Thus, the positive cone is positively invariant for 0 because, all 
trajectories that start in remain in R^ for all positive time. From system 0, with the initial 
conditions T(to) = Tq > 0 and G{to) = Gq > 0, we obtain the following system 

dG 


< (tg - ^go)G - fJ-cG'^, 

< (tt — dT)T — 

T{to) = To, 


dt 

dT 


G{to) = Go. 

Using the maximum principle, we deduce that 

G < 


Go 


^ + (l - ^ ) expi-Xofict} 


T < 


To 


n 

Yt 


+ 1 - 


7o 

Yt 


exp {—YTUrt} 


When t —>■ oo, we obtain 


lim G{t) < Xg = 

t^OO 


IG - < 5 go 


TG 


7G 


lim T{t) <Yt = 

t—>-oo 


ryrp — ^rp 

flT 


= AVI 1-7 

IT 


( 12 ) 


(13) 


(14) 


Hence, when 7 ^ > (ico 77 ^ > 6t, all trajectories of system 0 that reach the neighbourhood 
of B converge inside as t tends to infinity. Since B C R^, then B is positively invariant and 
attracting for system 0- 
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Appendix B: Proof of theorem 3.1 (existence of the semi-trivial pe¬ 
riodic equilibrium) 

Let T{t) = 0, from system Q, we have the following simple logistic impulsive differential system: 


^ = roG-ncG'^, 

, G(t+) = G{tn) - XfcGitn) t = tr. 


(15) 


The solution of system (15) is given in m- Here, is shown in detail the proof. 

Setting X = — in the first equation of system (15), we have the following differential equation 

Gr _ 


integrating 


from nr to t, we obtain 


dX 


dX 
-M = 


X(t) = ae-^G(t-nr)^ 

Using the variation of the constant a, we have the following differential equation 

da 


dt 




Thus, we have 


l{t) = ^e^ait-nr) ^ ^ ^ 

Xg 


Substituting © in (|16[), we obtain 


(16) 


(17) 


X{t) = — + b G R. 

Xg 


Considering X{t) = ~q^ (18), we have 


this implies that. 


At t = nr, we have 


^ 6 G R 

G{t) Xg ’ ’ 


Git) = 


1 


Xg 


b G R. 


G{nT) = 


1 


h G R, 


Xg^^ 


(18) 


(19) 
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and therefore 


1 


b = 


1 


G(nr) Xq' 


Substituting the expression of 6 in (19), we obtain 

G{nT)e^G{t-nT) 


G{t) = 


1 + 


G{nT) 


We have 


Xg 

G(f++i) = G((n + l)r) = 


^^rcit-nr) _ 

G{nT)e^°'^ 

-^G 


By substituting (21) into the second equation of (15), we obtain 

(1 - \fG)G{nT)e^Gr 


G{nT) = (1 - \fG)G{{n + l)r) = 


which implies that 


G{nT) = 


^G 

^G{(l-A/G)e^g"- 1 } 

e^GT _ 1 




Substituting the expression of G{nT) into (20), we have 


(20) 


( 21 ) 


( 22 ) 


G[t) = 


Xg{{1 - A/G)e’'G^ - i}e'^G{t-nT) 


{(1 - XjgV^^ - l}(e^G(i-nr) 


{e^GT - 1 ) 1 + 


(e'-G-r _ 1) 


Xg{(1 - XfGy^^ - Ije^Git-nr) 

{e^GT - 1) + {(1 - XfG)e^Gr _ l}(e'-G(t-nT) _ 1) 

Xg{( 1 - A/G)e’'G^ - Ije^Git-nr) 

{(1 - A/G)e'’G-^ - l}e^G(t-nT) XfGC^GT- 


Thus, when > 1, there exists 


Ge{t)=XG- 


{{1 - XfG)e^Gr ^^rait-nr) 


{(1 - XfG)e^GT _ i]e'^G{t-nT) _g A/gc^’g’ 

This completes the proof. 


> 0, f G [nr, (n + l)r[, n = 0,1, 2,.... (23) 


Appendix C: Proof of theorem 3.2 


We calculate the unique periodic equilibrium. 

Set X = — in the second equation of (1). We obtain the following differential equation 


dX 

dt 


—ttX + ht- 


(24) 
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Integrate system (24) in nr < t < (n + l)r, we have 


T{t) = 


1 + 


r(nr) 

Yt 


^rrit-nr) 


- 1 


Now, we solve the first equation of 0 

dG 


dt 


= ilG — ^go)G — iigG'^ — ■^tgTG, 


where the expression of T is given by (25). 

Set Y = — in (26). We have the following differential equation 

Gr 1_I 

dY 

— = -[re--fTGT{t))Y + fiG- 


Integrating 


from nr to t, we have 


dV 

_ = _^rG-7TGT{t))Y, 


ln(F) = —rG{t — nr) + 'Jtg / T{u)du + a, a G R 


= —rG{t — nr) + In 
d'T 


Tinr) 

1 + 

prrd-nr) _ .| 1 

Yt 

L -1 _ 


n, n G R, 


which implies that, 


where. 


Y{t) = Pe-^d,nT,T{.nT)) ^ PgR, 


m{t,nT,T{nT)) = rG{t — nr) + In 

Ht 


1 + 


r(nr) 

Yt 


p-Td-nr) _ ^ 


Variation of P gives 


where. 


Then, we have 


P{t) = dG fnrX{u,nT,T{nT))du + b, 6 gR 

x(Tnr,T(nr)) = ^ 

P{t) 


Y{t) = 


x{t,nT,T{nT))' 


(25) 

(26) 

(27) 

(28) 


(29) 

(30) 

(31) 

(32) 
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which implies that 


G{t) = 


From (25) and (33), we have 


' G{t) = 
T{t) = 


G{nT)x{t,nT,T{nT)) 

1 + ncGiriT) x{u, nr, T{nT))du 


G{nT)x{t,nT,T{nT)) 

1 + ncGinr) fl^x{u,nT,T{nT))du 


T{nT)e^T{t-nT) 

T{nT) _ 

Yt ^ 


1 ) 


(33) 


(34) 


where G{nT) and T{nT) are values of grasses and trees biomasses respectively, immediately after 
the pulse of fire at the time nr. G(nr) and T(nr) can be viewed as the initial values of Q in the 
interval [nr, (n + l)r[. The initial values may change in different intervals. For all t G [nr, (n + l)r[, 
using the fact that. 


G{tt) = G{tn)-\fGG{tn), 

^(^n) = T{tn) — XfTUj{XfGG{tn))T{tn), 
we have the following discrete system for G^nr) and T{nT) 

(1 - XfG)G{nT)x{{n + l)r, nr, r(nr)) 


G((n + l)r) = 


r((n + l)r) = 


Setting 


U{G{nT),T{nT)) = 


y(G(nr),T(nr)) = 


1 + HGGinr) x(n, nr, T{nT))du ’ 

(1 — Ajrti;(AjG'G(nr)))r(nr)e'’^’" 

(1 - A/G)G(nr)x((ra + 1 ) t , nr, r(nr)) 
1 + HGGinr) x{u,nT,T{nT))du 

(1 — A/T<^(G(nr)))r(nr)e^^’’ 


implies that system (36) is equivalent to 

G((n + l)r) = U{G{nT),T{nT)), 

r((n + l)r) = ■F(G(nr),r(nr)). 


(35) 


(36) 


(37) 


(38) 
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The existence of a periodic solution of Q (with period r) is equivalent to an existence of the 
equilibrium of the discrete system (38). This leads to solve the following system 


We have 


U{x,y) = X, 
V{x,y) = y. 


U{x,y) 


X 1 + ycx 


f(n+l)T 


x{u,nT, y)du 


(1 - \fG)x{{n + l)T, nr, y), 


which implies that 


(39) 


{1-XfG)x{(.n + l)T,nT,y)-1 
X = -r-T-r^- := (i)[y). 


HG xiu,nT, y)du 


(40) 


On the order hand, we also have 


yix,y) = y 4^ y = 


Yt{{ 1 — XfTUj{XfGx))e'^x'r _ i| 

{e^TT _ 1 ) ■ 


which implies that 


^(e-xr _ 1 ) + XfTU;{XfGx)e^xr ^ 

Yt 


(41) 


From (41) and the nonnegativity of the variable y and the function uj{XfGx) with x > 0, it 
follows that y must belong to the interval D = [0, Ir]- Substituting x = 4>{y) into the left side of 


(41) yields an equation for y 

y 


Yt 


(e-xr _ 1 ) + XfTOoiXfGmV^^ - - 1 ) = 0 . 


Set 


y 


h{y) = - 1 ) + XfTOj{XfG<t>{y)V^^ - - !)• 

Yt 


(42) 


It is obvious that h{y) is nonnegative and continuously differentiable with respect to y. The 
algebraic calculation shows that 


We have 


h{0) = XfTu;{XfGm)e^^^-1), 

hiYr) = XfTUj{XfG4>{YT))e'^x'r^ 

(1 - XfG)xi{n + 1)t, nr, 0) - 1 


(43) 


m = 


d-G xiu, nr, 0)du 


where 


'■(n+l)r 


x((n + l)r, nr, 0) = and 

»(n+l)T 


xiu,nT,0)du= / = — (e^G^-l). 

J nr 
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Then 


m = ^ = Geir). 


Note that 


eroT _ 1 
yGe 


Geir) >0 4^ '^0,pulse > 


From (431, it is easy to see that 


hiYr) > 0 . 


We also have 


(44) 

(45) 


Then 


h{0) <0^ XfTU}{XfG(l){0))e'^^^ - - 1 < 0 

^ 1 < - XfMXfGGeir))) 

c^T-r 




1 


1 — XfTUj{XfGGe{T)) 

1 


< e 


<t4> In ^ 
4 ^ 1 < 


1 - XfTUj{XfGGe{T)) ^ 

rr 


< TtT 


-In 

T 


1 


1 - XfTUj{XfGGe{T)) 


4^ 1 <TZ, 


0,pulse ‘ 


h{0) < 0 T^Q^pulse > 1- 


(46) 


Thus, when (45), and (46) are verified, there exists at least one positive zero of h{y) = 0 in the 
interval D. To have the uniqueness, we show that h is a highly monotone function. 

The derivative of h{y) with respect to y is 

^ - 1 ) + <p'{y)XfTXfGU^'iXfGcPiyW^\ 

where 


- </> 2 (y)> 

with 

s i'i-- >^fG)x {{n + l)T,nT,y)yG xiu,nT,y)du 

^- r(n+i)r , -TTT^-’ 

\yGjnr x{u,nT,y)duy 

and 

. ((1 - A/G)x((n + l)r, nr, y) - l)/iG x{u, nr, y)du 

^- Kn+I)r , -TTT^- 

yyGjnr x{u,nT,y)duy 
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We have 


X {u, TIT, y) = 


^^VTiu-nr) _ 


'iTG _ Yx _ m(w,nr,y) 

/iT 1 + _^(^^rT{u-nT) _ 

Yt 


ITG {e^T(u-nr) _ 

Tt \ _^t^rT{u-nT) _ 


—x{u,nT, y). 

1 ) 


Substituting x {u, nr, y) in (p'i{y) and 02(?/)) we obtain 


(1 - ^/g) < 


My) = 


ITG - 1) 

i + M^e'^TT _-^y 

Yt 


-X((n + l)r, ITT, y) 


where 


MG nT, y)du 

= -(1 - XfG)x{{n + l)r, riT, y)0, 

( XTG - 1 ) \ 

1 + A(e'-Tr_i) 
_£T_ 

yG x{u,nT, y)du 


e = 


\ 


and we also have 


-{(1 - XfG)x{{n + l)r, riT, y) - 1} 


Ttg 

Tt 


(n+l)T ^^rxiu—nr) _ 


1 ) 


My) = 


1 + A(e^T(«-nr) _ 

Yt 


x{u,nT, y)du 


(mg M: riT, y)du) x{u, riT, y)du) 

[e^Tiu-nr) _ M{^(TTr _ ; l )} 

-{(1 - XfG)x{{n + 1 )t, ut, y) - 1} Z -“T- v , S.. - —x{u, tit, y)du 


= e- 


{e^TT _ 1){1 + A(e"T(«-nr) _ 

Yt 


Setting 


Tr{u, UT, y) = 


Xiu, riT, y)du 


[e^Tiu-nr) _ _ ;l)} 

_£T_ 

{e^TT _ 1){1 + ]L{^yrT{u-nT) _ ^ 

Yt 


we obtain 


4>{y) = (t>i{y) - 4>2{y) 
= &Q{y), 
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where 


Q{y) = 


{(1 - A/g)x((« + l)^, nr, y) - 1} vr('u, nr, y)x{u, nr, y)du 


From the expression of 


X{u, nr, y)du 

dh{y) 


-{l-XfG)x{{n+l)T, nr, y). 


dy 


we see that that h{y) is monotonous for all r. Hence, there is only 


one positive root of h{y) = 0 for r > 0. Therefore, when T^Qp^ise > ^ pulse > equation 

h{y) = 0 has a unique positive zero for r > 0. Thus system Q admits a unique non trivial periodic 
solution. 


Appendix D: Proof of theorem 3.3 (Local stability of constant equi¬ 
libria). 


The proof of the stability is on the basis of the linearization to (38). Letting (G'^,r®) be the 


linearized system of (38) is 


where. 


with 


We have 


r (0;yT)- Set X{t) = G{t) — G®, and Y{t) 

X((n + l)r) 

= 

aiiX(nr) + ai2T(nr), 

y((n + l)r) 

= 

a2iX(nr) + a22T(nr), 

dU , , 


dU, , 


Te) 

, Ol2 = ^(G^n, 

dV , , 


dV, , 

»21 = ^(0 

T^) 

, 022 = ^(G^n, 


U{x,y) = 


V{x,y) = 


(1 - XfG)xx{{n + l)r, nr, y) 

1 + HGX x{u, nr, y)du ’ 

(1 - XfTOj{XfGx))ye^TT 

l^A(e-rr_i) 

Yt 


(1 - A/G)x((n + l)r, nr, T^) 

On — -TTn- ^ 

(l + /iGG®X!” ^^x(n,nr,re)dn)2 


021 = 


'~pe 


< 0 . 


(47) 


(48) 
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Setting ^ = 


— 1 ) 




< 0, we obtain 


012 = 


(1 - XfG)G^x{{n + l)r, nr, |l + /igG® (1 “ nr, T^)) x(n, nr, r®)(in| 

{1 + ^igG^ x{u,nT,T^)duY 


where, 


7r(n, nr, T®) = 


(grT(«-nr) _ —{e^Tr _ 

Yt 

(gr-rr _ 1){1 + ^(grr(«-nr) _ 

Yt 


It is easy to show that for nr < n < (n + l)r, we have 7r(n, nr, T®) < 1. Then, ai 2 < 0. 
We have, 

(1 — Xf'j'Ui^Xj'GG^})^^^^ 


022 = 


l + |(e--l) 


> 0. 


The stability of the equilibrium of (38) can be determined by eigenvalues of the linearized matrix 
A = (ajj) 2 x 2 - 

1. The two eigenvalues for A at (0; 0) are 

pi = (1 - XfG)x{{n + l)r, nr, 0) 

= (1 - XfG)e^^\ 

and 

P2 = > 1. 

Then, the trivial equilibrium (0, 0) is always unstable. 

2. Concerning the local stability of forest equilibrium, we have the following two eigenvalues 


ni = (1 - XfG)x{{n + l)r, nr, Yr) 
= (1 - A/G)e^ 


ITG 

tqt— - rj^T 

Pt 


= (1 - A/G)e 


rarl 1 - 


1 


T^oi 


where 


The second eigenvalue is 


7^ol = ^x 


rx 'jtg 


V2 = - < 1 . 

grrr 
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Thus, we have the following results: 

• If 7^01 < 1) then ui < 1, which implies that the forest equilibrium Eqi = (0;Tr) is locally 
asymptotically stable (LAS) (similarly in the continuous model (see Tchuinte et al. 2014 

m))- 

• If 7^01 > 1 and ■^o,7^ol = pulse “ '^) ^ is LAS since ui < 1. This 

situation is specific for the impulse model. The continuous model does not imply the stability 
of the forest equilibrium when T^oi > 1- 

• If 7^01 > 1 and 7^o,7^ol > (0; ^r) is unstable. 


Appendix E: Proof of theorem 3.4 (Local stability of the periodic 
grassland equilibrium) 


To show that (Ge(t), 0) is LAS, we consider Floquet’s theory. Set G{t) = Ge{t) + x{t), and T(t) = 
0 + y{t), where x{t), and y{t) are small perturbations, every solution of the linearized equations can 
be written as 


x{t) 

y{t) 


= m 


x(0) 

1 /( 0 ) 


where, 4>(t) = ,i,j = 1,2 is a fundamental matrix and satisfies 

d^{t) 


dt 


= A{t)^{t), 


(49) 


with. 


A{t) = 


XG - ‘^fJ-cGeit) -'iTcGeit) 


0 


XT 


Since 4>(t) is the principal fundamental matrix, then 4>(0) = I 2 , where I 2 is the identity matrix 
of 7W2(R). Then, 


x(nr+) 

y(nr+) 


1 - A 


/G 


0 


x(nr) 


0 1 - A/Ta;(A/GGe(nT)) / \ 2 /(nr) 


and hence, if the absolute value of all eigenvalues (Floquet multipliers) of the monodromy matrix 


M = 


1 - XfG 0 

0 1 — A/rw(A/G(5e(r)) 


4>(r) 


are less than one, the periodic solution (Ge(t);0) is locally asymptotically stable. 
By calculation, we obtain 
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m = 


0 (^22 (t) 


where, 


(^11 (t) = exp \ ret - 2hg / Ge{u)du ^ , 


and 


1 . 


(P 22 {t) = exp I J ttcIu^ 

We deduce that, eigenvalues A*, i = 1, 2 of M are 

Ai = (1 - A/G)</?ii(r), 

A 2 = (1 — \fT‘jj{^fGGe{T)))(p22{T). 

From the hrst equation of Q, we have 

= [tq- - ^TGT)dt. 

Cj 


Integrating (51) in [0,t], we obtain 


G{t) = G(0) exp jrar - I^G j G{u)du — jtg j T'(rt)drt| • 


At (Ge(t);0), equation (52) becomes 


Geir) = Ge{0) exp <j r^r ~ I^G J G{u)du 


From the third equation of (Q, we have 


Ge(0) = (l-A/G)Ge(r). 


Substituting (54) in (53), we obtain 


which implies that 


Ge(r) = (1 - \fG)Ge{T) exp tgt -^iG G{u)du 


(1 - A/G)exp <1 rcT - G{u)du\ = l. 


(50) 

(51) 

(52) 

(53) 

(54) 
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Thus 


Ai = (1 - A/G)<y7ll(T) 

= (1 - A/g) exp jrGT - 2/iG j G{u)di 


= (1 - A/g) exp < roT -HG G{u)du ^ exp -/xg / G{u)d' 


/o 


= exp —/xg J G{u)du f < 1- 


XXX 


On the other hand, we have 

A 2 < 1 Tlo^pulse 


rx 


< 1 . 


-In 


^ \1 — A/Tw(A/GGe(r))y 

Then, the periodic grassland equilibrium Eiq is LAS if ^ This completes the proof. 


Appendix F: Proof of theorem 3.5 (local stability of the periodic 
savanna equilibrium) 


Now, we investigate local properties of the periodic savanna equilibrium. Similarly to the proof of 
the local stability of the periodic grassland equilibrium, we set G{t) = G*{t) + x{t), and T(t) = 
T*{t) + y{t), where x{t), and y(t) are small perturbations and they are solutions of the linearized 
equations 


/ x{t) \ / x(0) \ 


(t) is the fundamental principal matrix which satisfies 


d^*{t) 

dt 


re -‘^IJ'gG* (t) -'JtgT* (t) -^TGG*{t) 

0 rx — 2/j,xT*{t) 


By calculation, we have 




where. 


^*{t) 


0 ¥^22 (^) 


(fliit) = exp ( / (rG - 2yLGG*{u) - 'yxGT*iu))du 
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and 


(^22(i) = expjy (rr - 2/XTr*(u))du| . 

From system Q, at t = tn, we have the following system 

" - )(■” 

y(nr+) / \ 0 1 - A/Ta;(A/GG*(nr)) / \ y(nr) 

Hence, according to the Floquet theory, if all eigenvalues (Floquet multipliers) A^ and A 2 of 

1 - Xfc 0 


M* = 


0 1 — A/tw(A/gG (r)) 


^*ir) 


are less than one, then the coexistence Tree-Grass periodic equilibrium is locally asymptotically 
stable. We have 

Ai = (1 - A/G)‘^n(r) = (1 - Xfc) exp jrcr - 2/iG £ G*{u)du - ^tg fo f*(u)duj , 

\* = {1 - XfTUj(XfGG*(T))}<fi^ 2 M = (1 - XfTUj(XfGG*(T)))exp jrrr - 2 ^t/J" f*{s)ds^ . 

(55) 

Starting with A^, we have 


A^ < 1 (1 — A/g) exp tgt — 2^g I G*{u)du — jtg J T*{u)du j < ^ 

-In ^ ^ 


^1_ 2/xg 1^1 _jTGfl r f*^u)dv\ < — ^ 


rG \r Jo 


fG \T Jo 


vg 


r G*{u)du] - ^ rf*{u)du]< ^ 


Xg Vr 


^G \T Jo 


n: 


Ge 

O^pulse 


Integrating G*{t) and T*{t) from 0 to r, we obtain 


f^G*{u)du = — IuIi + ugG* / xiu,0,T*)du] , 


d-G 


£f*iu)du = ^ln(l + ^(e'--"-l) 
fJ-T V 


Substituting (|^ which is expression of T* in the second equation of (56), we have 


[ T{u)du = — [ln(l — XfTOj{XfGG*)) -|- r^r] 
/o d-T 


(56) 


(57) 


42 











Then, 


Ttg f 1 


re \T Jo 


f*{u)du] = 


7TG 

1 

VG 

d'T 

ITG 

rx 

VG 


1 


'Tim 



vt -In 


1 - 


T {1 - XfTU}{XfGG*)) 
1^^ 1 
T (1 - XfTU}{XfGG*)) 
VT 


1 - 


7^; 


:=7^: 


0,stable , 


G* 

0,7^01' 


Thus we have 


2 r \ ^G* 1 


< 1 ^ 1 - ^ / G*{u)du ) < + 


Ge 


. ^G* 1 2 

^ 1 < ^o.Teoi + 


K 


O^pulse 


1 


V' 

^O^pulse 

G*iu)du) ■.= n*o%able- 


Therefore, 

xi<i^n*o%,,i,>i. 

We have x(7‘, t", T*) = 1, which implies that 

x{t,t,T*)G* 


(58) 


Then, 


G*(t) = 


1 + /xgG* x{u, r, T*)du 
= x{r,T,T*)G* 

= G*. 

A 2 = (1 — XfTbo[G*)) exp {—rj^T — 2 ln(l — Aj 7 ’Cj(AjgG*))} 
= exp{—r-pr — ln(l — XfTUtiXfGG*))} 

= exp I —r^T + In 


= exp < 


/i 


VtT 


In 


(1 — XfT<jj{XfGG*)) j 

1 \ 


r (1 - XfTUj{XfGG*)) 


- 1 


V 


rx 


= exp < rxT 


7^; 


- 1 


0,stable 


Thus, 


A 2 < 1 


'^0,stable > 


(59) 


From (58) and (59), it follows that the coexistence Tree-Grass periodic equilibrium is locally 
asymptotically stable if ”^0 siaftie ^ ^ '^*0 stable ^ This proof completes this section. 
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(Global stability of the forest 

equilibrium) 

From system Q, we have 


Appendix G: Proof of theorem 3.6 


G'{t) < rcGit), t^tk, 

Gi4) = (1-A/g)G(4) t = tk. 

Using impulsive differential inequations (Lakshmikantham et al 1989 |58|1. we have 


(60) 


G(t) < I (1 - A/g) ) exp r / reds) 


' t 

= ^ 0(1 - A/g) '^1-’’' 
< Go(l — A/g) 


io' 

^ ^ exp{rG(t - to)} 

W A 

^-1/ exp{rG([t - to] + 1)1 


t 


= Goe^‘='^(l - A/g) AU . 

(\* 

= Goe^^'^{l - A/g) 


4 

^ ^ exp{rG([t - to])} 
to' A 

/ exp < rcT 


T 


= Goe^cr g^p , 


rcT 


1 - 


-In (-— 

T yi-A/G 

re 


Y 

A 


= exp < tgt 1 — 


n: 


0,pulse ^ 


Then, T^^4,ise ^ ^ implies that G(t) —)■ 0, as t —)■ 00 . 
Now, we will prove that 


lim r(t) = Yt. 


AVhen G(t) = 0 and nr < t < (n + 1 )t, n = 1, 2, ...,N-r, system 0 becomes 


= rrT - htT'^, t^tn, 

dt 

. T{t+) =T{tn) - XfTl^{0)T{tn), t = t, 


(61) 


(62) 
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(63) 


since cu(0) = 0, system (62) is equivalent to 


dT ^ ^2 

— = rrr - ,.Tr . 

which has two equilibria 0, and Yj’. It is obvious that is GAS. Thus, 


lim T(t) = Yt. 

t^+oo 


This complete the proof. Thus, when T^Qp^ise < ^ forest equilibrium Eqi is GAS. 


Appendix H: Proof of theorem 3.7 


We prove the global stability of Eio in the following two steps: 


Step 1. First, we show that lim T(t) =0 if 


t^OO 


Ge 

0,stable 


< 1. In fact, from system (11, we obtain 


(64) 


T'{t) < rTT{t), 

= (1 — A/Tw(A/GG'e(tn)))T(tn), t = tn- 

Using impulsive differential inequations (Lakshmikantham et al 1989 |58]) we show that 

T{t)<T{t'^)( {1 - XfTUj{XfcGe{tn)))\e^p ( f rxds] . 

\to<tn<t j J 

It is obvious that Ge increases monotonically in [nr, (n + l)r[, n = 1, 2,... On the other hand 
a; is an increasing function, then for all n = 1, 2,... we have 


(1 — XfTOj{XfGGe{nT))) < (1 — A/Tw(A/GG'e(r))). 
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Then, we have 


T{t) < T{t^) ( n “ -^/Tw(A/GGe(T))) j exp ( f rrds^ 


Kto'^tk 

= 7o(l — \fT0j{\fGGe{T))) 

< To(l - XfT0j{XfGGe{T))) 

= Toe^-^{l-XfTOj{XfGGe{r))) 

= Toe^^ exp < r^r — In • - 


■ t' 


to 

_T 


_ T _ 

■ t' 


to 

T 


T 


exp{rT(t - io)} 
exp{rT([t - to] + 1)} 


exp < ttt 


1 


= Toe'"^ exp < 


iln 


TtT 


1 — XfTto{XfGGe{T)) ^ 
1 


1 - 


^ \1 — A/rw(A/GGe('J‘))^ 

rr 


\ 


/ 


= exp < vtt 1 — 




0,pulse 


Thus, when TZ, 


0,pulse < T{t) ^ 0 as t ^ oo. 


Step 2. We prove that lim |G(t) — Ge{t)\ = 0. 

t^OO 

Since lim T(t) = 0, then for any given ei > 0, there exists ti > 0, such that 

t^OO 

-ei < T{t) < ei. 

For all t > ti. Using (65) into the first equation of ([^, we obtain 


dC 

— > tgG - hgG^ - yrceiG, t 7^ U 
at 

, G(t+) = G(t„) - A/GG(t„), t = tn. 


(65) 


( 66 ) 
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Let z = G we have — = Then system (|66|) changes into following 


dz 


1 dG 


dt 


G2 dt 


dz 

-^ < -[re - jTG^ljZ + HG, 


< ^ (67) 

^(tn) = -. -- ^(tn), t = tn. 

1 - A/G 

Set = tq — jtg^i- Using impulsive differential ineqnations (Lakshmikantham et al 1989 
ESI), we have 


’ ..E, (Evs) (" I 


= EE) 


l-A 


/G 




+ fJ-Gexp{-Aeit) 


Hi 


n 

S<tk<t 


exp - tl)) 


exp (^ejs) ds, 


l-A 


/G 


which implies that 
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z{t) < z{tt) 


1 - A 


fG 


h 


exp -ti)) 


+ ^Gexp(-A,^t) [ / TT ---— exp(A,js)<is 


'‘[^i]+2t 


'^[tl]+T S<tk<t 

+...+/'” n 


^ exp (ylgj s) ds 


'[t] 


1 - ^/g/ 


S<tj^<t 


= z{tt) 


1 - A 


/G 


+ -^exp(-yl,,t) [ 


[ii]+T 


h 


exp (-ylei(t - h)) 

n ( i - A/g ) d{A,,s) 


S<tj^<t 


+ 


|•[t-l_] + 2T 
[* 1 ]+"^ S<tk<t 

rlt] 


+ ... + 


'[t]-T 


n ( 1 —\—1 ®^p 

'l<i VI - a/g; 

S<tk<t 


+ 


'[t] 


1 - A/g/ 
n ( i - A/g ) 


S<tk<t 


= z{tt) 


1 - A 


/G 


^ei 


+ 


1 - A 


1 - A 


/G 


fG 

h' 


t — ti 






(^g^ei([G]+2r) _ gA£^([tl]+T)^ 


+ ... + f(e^'i - l)e^n(W-^) 
VI - A/g/ 


\/G 

+ e^'i*-e^^iW]. 
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Using the fact that [U] <ti, it follows that 


z{t) < z{tj) 


[rlJ 


t — ti 


1 - ^/G/ 
ti rti] ^ 
r r ) 




(^g^ei([*l] + 2T) _ g^ei([il]+T)^ 


+ ... + ( ^ ) (e^^i - 

Vf - A/g/ 

gA^ji _ gA,Jt]] 




I'Tjj (-A,Ut-ti)) 


+ M(eA.i^_l)e(-Ani) 


1-A, 


[((l-A^G)e^^iO'^ 


ff-1- 

+ ... + ((l-A/G)e^'i")^L^^ y 

_)_ i^^Lg(-Aeii)jgAejt _ g.4,Ut]j 


We have 


= (l-A/G)e^nU 


5 = aL^J +aVL^ 


+ ... + a 


= aL’^-l {l + a + a^ + ... + a”} 


= X 


1 - 


In the expression of S, we have — + n = — — 1, which implies that 

[rj [rj 

't] [ti] 
n = — — — — 1. 
r r 


49 



Substituting (69) in the expression of S leads to 


S = al 


pr 


't' 




'ti 


' t' 

r 

1 — a 
X - 

T ^ 


T 

a 

T _ 

— a 

(r J 


1 — a 


1 — a 


( 70 ) 


Using (68) and (70), we have 


z{t) < L{t) + < 


1 - (1 - A/G)e^’=U 


_|_ _ e ^'^1 W 


where 




\/G 






1-A 


fG 


(e^n-) 


1-(1-A/G)e^'i^ 


Then we have for all t > ti, 


z{t) < L{t) + < 

^ei 


1 - 


- 1 ) 


Ae, r 


1 - (1 - XfG)e^^i 




Since ei > 0 is arbitrary, it is obvious that 


lim L(t) = lim z(tt) -- 

t—>-00 i—^oo \ 1 — A 


Vg 


= (-Aei(i-ti)) 


+ lim 
i^oo A, 


1-A 


fG 


'h 

{e^nr)[T 


l-(l-A/G)e^^G 


= 0 . 
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Thus, it follows that 


( 71 ) 


lim z{t) < 

t-^oo ^ ~ 

On the other hand lim = lim ire — ^TG^i) = ^G- Thus for ei 

ei^o ei—^-O 

lim zit) < . 

*—>■00 Xq 


0, (71) changes to 


(72) 


Because z = —, from (72) it follows that 

Gr 


G{t) - Xg ’ 


which implies that 


lim G{t) > Xg- 

t^OO 


Thus for any €2 > 0, there exists t2 > 0 such that 


G{t) > Xg — £2 


for all t > t 2 - 

On the other hand into B, we have 


G{t) < Xg, 


for all t > 0. Coupling (74) and (75), for all t > t 2 'we obtain 


(73) 


(74) 


(75) 


> G{t) > Xg - 62 


Let 6 = min{ei, 62}, and 


max{ti,t2}- Then from (76), we obtain 


Xg > G{t) > Xg - 6, 


(76) 


for all t > tig that is 


It is obvious that 


0(t) — y Xg, ^rs t —y 00 , c —y 0”^. 

lim Ge{t) = Xg- 

t^OO 


Using (77) and (78), we have 


(77) 

(78) 


lim \G{t) — Ge{t)\ = 0. 

Therefore, the grassland savanna periodic equilibrium is globally asymptotically stable when 
'R-q pulse < 1- This completes the proof of the theorem. 
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